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Q I We show that in generic leptoquark (LQ) extensions of the standard model lepton 

and baryon numbers are broken at the level of renormalizable operators. In particular, 
this may cause fast proton decay unless the leptoquarks are heavy enough. We derive 

CN I stringent bounds for the 1st generation LQ masses and couplings from the proton stability 

' constraints. 
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O ! Leptoquark models [|I|, 0] remain an attractive possibility for new physics beyond the Stan- 

^ I dard Model (SM) admitting non-SM particles, leptoquarks (LQ), with masses at the elec- 
O ' troweak scale. LQs are vector or scalar particles carrying both lepton and baryon numbers 
^ ■ and, therefore, having a distinct experimental signature. For these reasons searching for LQs 

is a promising subject for present and near future experiments 0, ffl. 
QhI Theoretical motivation for LQ models usually refers to the low energy limit of some more 

^ I fundamental theory associated with an energy scale much higher than the electroweak scale. 
■ In the literature it has been argued that the origin of LQ models may reside in grand unified 
theories (GUT) p, [0], models of extended technicolour composite models |lOl and some 



X 



other high-energy scale model. However the arguments within this framework in favor of light 
5^ I LQs with masses of the order of the electroweak scale are quite vague. Nevertheless the generic 
structure of LQ models can be completely determined by the symmetry with respect to the SM 
gauge group p|, |ll|. In this way LQ models with light LQs can be studied without referring to 
their high energy scale origin. 

Adopting this approach we examine the question of Lepton (L) and Baryon (B) number 
conservation in LQ models. We show that in the LQ Lagrangian there exist renormalizable 
AL = 1 and AB = 1 terms which may affect the conventional LQ phenomenology. In particular, 
the previously overlooked in the literature AB = 1 terms may cause fast proton decay unless 
the LQs are heavy enough. Thus the proton stability constraint casts lower limits on the LQ 
masses. We show that these limits are more stringent than those existing in the literature [0, ^ . 
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Here we construct the interaction Lagrangian of a generic leptoquark model retaining all 
the renormalizable couplings invariant under the SM gauge group SU{3)c (8) SU{2)l ® f/(l)y. 
We separate the LQ interaction Lagrangian in the following three parts 

C-LQ-int = ^LQ-l-q + ^LQ-H + ^LQ, (1) 

corresponding to LQ-lepton-quark, LQ-Higgs interactions and LQ self-interactions. 

The LQ-lepton-quark interaction terms jCiq^i^q, which fix the LQ field content, were con- 
structed in Ref. for the scalar S and vector leptoquarks. Here we show only the scalar 
LQ interactions 

Cs-i-, = -IT-PRe-S^^ + \f^ ■WPRe-~Sl + \[f^-uPLl-S^]^ + \y2-d^^^ 

+ AS,'') • -^-PlIT^I ■ Sj^^ + ■ qPRiT^e ■ S^^ + Ai • q-^PLtr^Sll + h.c, (2) 

where Pl,r = (1 =F 75)/2; q and / are the quark and lepton doublets; S"] are the scalar LQs 
with weak isospin j=0, 1/2, 1, coupled to left-handed {i = L) or right-handed {i = R) quarks 
respectively (for a discussion on chiral couplings see Ref. [0). The LQ quantum numbers are 
listed in Table 1. For LQ triplets $i = Si, Vi we use the notation $i = r ■ $i. 

To obey the stringent constraints from FCNC processes it is usually assumed that the LQ 
couplings are generation "diagonal" , i.e. they couple only to a single generation of leptons and 
quarks. This implies that there exist three generations of LQs with the assignments of Table 
1. In general these couplings could involve all fermion and LQ generations. 

The renormalizable LQ-Higgs interaction terms for the scalar and vector LQs have been 
constructed in Ref. [|ll|]. Again, for brevity we show only the scalar LQ interactions: 

Clq-h = hfmT2~Si/2 ■ SI + h^HiT^S^ ■ Si/2 + Y^^ (mr^Sl/,) ■ {sI,h) + (3) 
+ Y, (Hit2S\h) ■ So + K^''> {h^S.h) ■ Si^ - [Ml - g^^'^ h) <l>*^t$^2_ 

Here H is the SM 5'f/(2)i-doublet Higgs field. is a cumulative notation for all the leptoquark 
fields with i = L,R (the same for ii,2)- 

The new, previously overlooked part of the LQ interaction Lagrangian, corresponds to the 
LQ self-interaction terms. For brevity we write down this part in the form of SM group singlet 
products of the LQ representations 

J^LQ = f^liSl^ X ^1/2 X + ^{Sl X Si/2 X 5i/2) + fi,{Sl X ^1/2 X ^1/2) + (4) 

+ gi,D ■ (K^ X Vi/2 X V1/2) + g'oD ■ {V^ x Vl^, x V1/2) + giD ■ {V,^ x V1/2 x VI/2) + 
+ r7,^($a X • ($^ X <l>^) + r]i\<i>}, X $1, X ${/2 x $1/2) + vl'iH x <l|/2 x $1/2 X 
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where ^ = S,V and the superscripts i,j = L,R are LQ chirahty indexes defined in Eq. (^. 
The subscripts a,P in the 1st quartic term denote all the types of vector and scalar LQs. 
The parameters /z^ and ga,Va dimensionful and dimensionless parameters of the model 
respectively. The SM gauge group covariant derivative D^, defined in the standard way, acts 
on all the LQ fields V in brackets. Here we suppressed the Lorentz indexes of the vector 
LQs which must be contracted in all the possible pairs to form Lorentz scalar products. Thus 
the above indicated vector LQ couplings represent groups of terms with different coupling 
constants for each term. For example, the first group of vector LQ couplings contains the term 
X Vi/2 X iDi,Vi/2) which has its own coupling constant fi^^^ etc. 

Let us check the lepton and baryon number properties of the terms in the LQ Lagrangian. All 
the LQ-lepton-quark interaction terms CiQ^i^q conserve both L and B numbers by construction 
and, therefore, are AL = /S.B = operators. In the sector of the LQ-Higgs interactions Clq-h 
the quartic terms are also AL = AB = operators while the trilinear terms violate lepton 
number but conserve baryon number so that in this case AL = 2, AB = 0. In the LQ self- 
interaction part Clq again the quartic terms are AL = AB = operators while the trilinear 
terms violate both lepton and baryon number as AL = —AB = 1. 

The trilinear AL = 2, AB = terms in the LQ-Higgs sector contribute to lepton number 
violating processes, such as neutrinoless double beta decay {OuPP) |jTl|, as well as to the Majo- 
rana neutrino mass matrix. Thus from the experimental limits on the rates of these processes 
and neutrino masses one can deduce the corresponding limits on the LQ model parameters. 

The trilinear AL = —AB = 1 LQ self-interaction terms in Eq. (^) can induce proton decay 
and, therefore, the existing stringent constraints on proton stability may produce valuable lower 
bounds on the LQ masses. These terms, combined with the LQ-lepton-quark couplings from 
Eq. (^, generate at low energies AL = —AB = 1 lepton-quark contact effective operators. For 
example, the 1st LQ self-interaction term in Eq. (^) generate the following dim=9 operator 

O(^) = XI III, ■ PRd){e Pnd) {u^ Pne), (5) 

where Ms is the typical mass of scalar LQs. This and the other possible AL = —AB = 1 terms 
can induce proton decay in the channels: 

p — > e~e~^u7r^7c~^ , e~uu7r~^7r^ , e^e'^vn^, vuvn^ . (6) 

The corresponding proton lifetime Tp can be estimated in the usual way. In this estimate 
we assume that all the dimensionless coupling constants in Eq. are of the same order of 
magnitude ~ A and that all the trilinear scalar LQ self-interaction terms in Eq. have 
the same mass scale /Xg ~ /ig ~ /^i ~ yU- We also assume that there is no strong cancellation 
between different contributions to the proton lifetime. With this assumptions we obtain for the 
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contributions of the scalar LQ the following estimate 

Here k is the dimensionless phase space factor. The typical energy scale of proton decay is the 
energy released in this reaction. In the above formula the energy scale is set by the proton mass 
trip. The difference between rup and the actual energy released in a specific channel of proton 
decay is absorbed in the factor k. From the existing lower experimental bound on the proton 
life time (channel independent) t^^^ > 1.6 x 10^^ years [jl2[ we obtain for the scalar LQ mass 



Mc> «:i/i2 /^^^V^' • 10 TeV ^ /^^V^' • 10 TeV. 



• iu lev ~ A ' 

VlGeVy VlGeVy 



Here we put k^/^^ ^ 1 which is a good approximation taking into account the very small 
exponent of the phase space factor. For the same reason this result is weakly dependent on the 
simplifying assumptions made before Eq. (^. In order to reduce the number of free parameters 
in the above formula and deduce more information on the lower bound for the LQ mass one 
needs some assumptions on the possible values of the coupling constants A and on the mass scale 
IX of the trilinear operators in Eq. (^). If the theory is in the perturbative regime one may assume 
that the dimensionless couplings A obey the condition A^/(47r) < 1. Assuming further that the 
LQ model originates from some GUT scenario one may also think that these couplings are of 
the order of the Higgs-fermion Yukawa couplings of the corresponding generation. However 
all the assumptions of this type crucially depend on the high energy origin of the LQ models. 
Therefore, following the common practice we keep in our constraints both A and Ms as free 
parameters. 

As to the mass scale /i of the trilinear operators in Eq. (^), it seems reasonable to assume 
that /i > Ai? ~ 250 GeV, where A^? is electroweak scale. This could be motivated by the 
observation that these operators are associated with physics beyond the SM whose typical 
mass scale is expected to be larger than the electroweak scale Kp. Thus, taking /x ~ Ai?, we 
obtain a "conservative" lower bound 

Ms > A^/^ • 25 TeV. (9) 

However, the actual value of the scale /i can be much larger than Ap easily reaching, for instance, 
the grand unification scale Mqut ~ 10^^ GeV. The latter case results in the constraint 

Ms > A^/2 • (5 X 10^) TeV. (10) 

Thus the typical constraint for the case of scalar LQs ranges between those in Eq. (|^) and in 
Eq. (0) 

Ms > A^/2 • (25 ^ 5 X 10^) TeV. (11) 
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The constraints for the case of vector LQs can be obtained directly from the Eq. (j^) by 
the substitution fi ^ g ■ rup, assuming that all the dimensionless coupling constants in Eq. 
(|^) are of the same order of magnitude Qq ^ Qq ^ Qi ^ g- This substitution is motivated by 
the observation that the energy scale of the trilinear vector LQ operators in Eq. (^) can be 
estimated as gx[energy scale of the derivative] and that the energy scale of the derivative is 
given by the mean momentum flowing in the LQ propagators. The latter is comparable with 
the proton mass rUp. Since the energy scale appears in Eq. with the small exponent 1/6, 
the difference between rUp and the actual energy scale is not important. Thus for the vector 
LQ mass we obtain 

Mv > X^^^g^'^ ■ 10 TeV. (12) 

There exist in the literature constraints on LQ models from accelerator and non-accelerator 
experiments (for a summary see, for instance, ^ IH])- The most stringent constraints for 

the 1st generation LQs follow from the measurements of Atomic Parity Violation and from the 
universality in leptonic vr-decays. The best constraints from these experiments are |^ 

Ms>\ - 3.5 TeV, My > \ ■ 6.6 TeV. (13) 

The comparison of these constraints with those in Eqs. (|D-(|lT]) leads us to the conclusion that 
the proton decay constraints for the 1st generation scalar LQs are significantly more stringent 
than other existing constraints. For the case of vector LQs the comparison of the proton decay 
constraints with the existing ones can not be made in a direct way due to the presence in 
Eq. (|12D of the vector LQ self-interaction coupling constant g. Nevertheless it is instructive to 
consider some sample values of the coupling constants A and g. For instance, if A ~ 1 the proton 
decay constraints comparable with the constraint in Eq. (|T3|) occur only for large values of the 
vector LQ self-interaction coupling constant > 0.1 which is unlikely. However for smaller 
values of the LQ-quark-lepton coupling A comparable proton decay constraints can occur at 
very small values of g. For instance, if A ~ 0.01 the corresponding value of this coupling is 
g ~ 10-^ 

In conclusion, we derived new constraints on the 1st generation LQ masses and couplings 
from the experimental lower bound for the proton life time. We have shown that for the case of 
the scalar LQs these constraints are more stringent than the corresponding constraints obtained 
from other experiments. 
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Table 1: The Standard model assignments as well as lepton L and baryon B numbers of the 
scalar S and vector leptoquarks (LQ). (Y — 2{Qf,^ — T3)). 



LQ 


SU{3), 


SUi2)L 


Y 


Q em 


L 


B 


^0 


3 


1 


-2/3 


-1/3 


1 


1/3 


^0 


3 


1 


-8/3 


-4/3 


1 


1/3 


'S'1/2 


3* 


2 


-7/3 


(-2/3, -5/3) 


1 


-1/3 


5*1/2 


3* 


2 


-1/3 


(1/3, -2/3) 


1 


-1/3 


Si 


3 


3 


-2/3 


(2/3, -1/3,-4/3) 


1 


1/3 


Vo 


3 


1 


-1/3 


-2/3 


1 


1/3 


Vo 


3* 


1 


-10/3 


-5/3 


1 


1/3 




3 


2 


-5/3 


(-1/3, -4/3) 


1 


-1/3 


V^l/2 


3 


2 


1/3 


(2/3, -1/3) 


1 


-1/3 




3* 


3 


-4/3 


(1/3, -2/3,-5/3) 


1 


1/3 
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